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ABSTRACT 

Head  injuries  cause  many  of  the  fatahties  produced  by  blast  effects.  In  such  injury  predictions,  the 
skuU  is  more  vulnerable  to  dynamic  pressure  than  to  overpressure.  Dynamic  pressure,  resulting 
from  the  blast  wave  sweeping  over  the  human  body,  causes  the  body  to  be  swept  along  behind  the 
wave  at  some  displacement  velocity.  Injuries  occur  when  the  moving  body  encounters  stationary, 
sohd  structures.  The  displacement  velocities  for  various  skuU  injury  levels  are  known.  Solving 
the  motion  equation  for  a  body  in  a  transient  flow  field  yields  the  displacement  velocity  as  a 
function  of  peak  overpressure  and  impulse.  Baker,  et  al.  and  Mercx  previously  attempted  this 
computation,  but  in  each  case  significant  deficiencies  occurred. 

In  Baker's  analysis,  both  diffraction  pressure  and  drag  loading  were  computed,  but  the  drag  on  the 
human  body  contained  no  adjustment  for  the  object's  motion.  Thus,  some  of  Baker's  displacement 
velocities  exceeded  the  particle  velocity  behind  the  blast  wave,  (a  physically  impossibihty), 
resulting  in  very  inaccurate  curves  at  low  peak  overpressures  (<  .4  psi). 

In  Mercx's  analysis  diffraction  pressure  loading  was  neglected,  and  in  calculating  the  particle 
velocity  behind  the  blast  wave,  Mercx  used  the  ambient  air  density  instead  of  the  density  behind 
the  wave.  These  two  deficiencies  cause  considerable  inaccuracy  at  overpressures  above  ~10  psi. 

In  the  current  Improved  Displacement  Velocity  (IDV)  analysis,  the  effects  of  both  diffraction 
pressure-loading  and  drag-loading  were  considered,  with  allowance  for  the  body's  displacement 
velocity,  and  with  the  correct  air  density  used  in  the  particle  velocity  equation. 

Pressure-impulse  diagrams  were  calculated  based  on  all  three  methods  for  four  displacement 
velocities  (corresponding  to  four  fracture  probabihties).  At  low  pressures  the  IDV  curves  closely 
match  Mercx's  while  lying  well  above  Baker's.  Between  3  and  10  psi  the  IDV  curves  generally 
agree  with  Baker's,  as  Mercx's  drift  to  the  right.  The  IDV  curves  extend  to  pressures  as  high  at 
1000  psi,  while  Baker's  curves  end  at  ~10  psi,  and  Mercx's  at  -100  psi.  A  comparison  of  results 
has  revealed  that  the  IDV  model  represents  a  bridge  between  the  two  earher  models,  eliminating 
the  major  deficiencies  associated  with  each. 


Report  Documentation  Page 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

AUG  1998 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-1998  to  00-00-1998 


4.  TITLE  AND  SUBTITLE 

Comparison  of  Head  Injuries  as  Predicted  by  Three  Different 
Displacement  Velocity  Models 

6.  AUTHOR(S) 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Engineering  Analysis  Inc, 715  Arcadia  Circle, Huntsville, AL, 35801 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 

See  also  ADM001002.  Proceedings  of  the  Twenty-Eighth  DoD  Explosives  Safety  Seminar  Held  in  Orlando, 
EL  on  18-20  August  1998. 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OE  PAGES 

19a.  NAME  OE 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


1.0  INTRODUCTION 


Most  analyses  dealing  with  the  blast  effects  produced  by  explosions  are  primarily  concerned 
with  the  prediction  of  damage  to  inanimate  structures.  In  many  cases,  however,  the  ultimate 
purpose  of  such  analyses  has  been  the  prevention  of  injury  to  personnel.  For  this  reason,  as 
part  of  a  company- funded  research  program,  a  hterature  survey  was  conducted,  dealing  with 
the  vulnerabUity  of  the  human  body  to  blast  [1-12].  As  part  of  this  research  special  attention 
was  given  to  the  human  skuU.  Two  models  for  predicting  blast  injuries  to  the  skuU  were 
identified  [7,  10].  Comparison  of  the  predicted  injury  levels  produced  by  the  two  models 
revealed  significant  differences.  A  careful  review  of  the  two  analytical  approaches,  including 
contact  with  the  authors,  or  co-authors,  confirmed  certain  deficiencies  in  each  technique.  To 
resolve  the  differences  in  the  two  previous  models,  the  Improved  Displacement  Velocity  (IDV) 
model  was  developed  by  EAT 

2.0  DEVELOPMENT  OF  IMPROVED  DISPLACEMENT  VELOCITY  MODEL 


In  calculating  personnel  injury  due  to  blast,  the  vulnerabihty  of  the  skuU  is  of  special 
importance.  As  the  blast  wave  envelopes  the  human  body,  the  body  begins  to  be  swept  along 
at  some  displacement  velocity.  Under  such  conditions  injury  to  the  skuU,  Uke  injury  to  other 
bones,  is  considered  to  result  primarily  from  the  impact  with  a  soUd,  stationary  object.  For 
this  reason,  skuU  vulnerabiUty  is  generaUy  expressed  as  a  function  of  displacement  velocity,  as 
shown  in  Table  1. 


Table  1.  SkuU  VulnerabiUty  [3,  4,  6] 


Displacement 
Velocity  (ft/sec) 
10 
13 
18 
23 


Fracture 
ProbabiUty  (%) 
0 
1 

50 

99 


Calculation  of  the  maximum  displacement  velocity  is  thus  the  fundamental  step  in 
establishing  skull  vulnerahility.  The  basic  problem  can  be  described  as  the  calculation  of  the 
motion  of  an  immersed  object  in  a  transient,  compressible,  one-dimensional  flow  field.  The 
basic  governing  equation  can  be  written  as 

m  dV/dt  =  -I-  (1) 

where 


m  =  mass  of  human  body 
V  =  displacement  velocity  of  human  body 


t  =  time 


Fdrag  =  drag  load 
Fdiff  =  diffraction  load 


Now, 


Fdrag  +  Fdiff  -  Pr  ^ 

(t<ts) 

(2) 

Fdrag  =  CD  Ps  (U^  -  V)  1  -  V  1  A  /  2 

(t>ts) 

(3) 

Ks  A 

(ts  <  t  <  td) 

Fdiff  =  Si 

(4) 

lo 

(t  >  td) 

where 


Pj-  =  reflected  overpressure 
Cp)  =  drag  coefficient  (1.3)  [13] 


Ps  =  gas  density  behind  shock 


Ug  =  gas  velocity  behind  shock 


Pg  =  gas  overpressure  behind  shock 
A  =  projected  area  of  human  body 
tg  =  rise  time  to  reach  stagnation  pressure 
t(j  =  rise  time  to  reach  peak  overpressure 
According  to  compressible  fluid  flow  fundamentals  [14], 


Pr  =  2Ps  (7  Poo  +  4Ps)/(7  Poo  +  Pg) 

(5) 

Ps  ~  Poo(2  Pqo  +  6  Pg)  /  (7  Poo  +  Pg) 

(6) 

Ug  =  5  Coo  Ps  /  Poo  (7  Poo  +  6  Pg) 

(7) 

Coo  —  Too 

(8) 

where 

Pqo  =  ambient  pressure 
Coo  =  ambient  speed  of  sound 
y  =  specific  heat  ratio  for  air  (1.4) 
R  =  specific  gas  constant  for  air 
Too  =  ambient  air  temperature 
According  to  Glas  stone  [15], 


^s  ~  ^o  ( ^  t  /  tp) 


(9) 


tg  =  3  w  /  2  /  Ug 

(10) 

td  =  (b  -1-  4w  /  2)  /  Ug 

(11) 

where 


Pq  =  peak  overpressure 
tp  =  positive  phase  duration 
w  =  width  of  human  body 
b  =  thickness  of  human  body 

Now,  if  Pg  varies  hnearly  with  time,  the  duration  of  the  positive  phase  can  be  expressed  as 

tp  =  ^Ipp/Po 

where 

Ipp  =  impulse  for  positive  phase 

Inspection  of  the  preceding  equations  reveals  that  for  constant  values  of  m  (mass),  A 
(projected  area),  w  (width),  and  b  (thickness)  of  the  human  body  the  displacement  velocity,  V, 
can  be  expressed  as  a  function  of  the  two  variables,  Pq  and  Ipp.  Thus,  the  displacement 
velocity  can  be  expressed  in  the  form  of  a  pressure-impulse  (P-I)  diagram. 

As  part  of  two  earher  studies  [7,  10],  attempts  were  made  to  solve  the  problem  of  calculating 
the  displacement  velocity  of  the  human  body.  Because  of  significant  discrepancies  in  each 
attempt,  however,  the  results  are  not  in  agreement.  In  the  first  attempt.  Baker  et  al  [7],  in 
calculating  drag,  did  not  take  into  account  the  changing  displacement  velocity  as  a  function  of 
time.  Baker's  model  is  especially  important  because  it  has  been  accepted  by  the  Department  of 
Energy  [12]  and  is  commonly  used  by  the  Department  of  Defense  [16]. 


Mercx  [10]  used  an  unreahstically  low  drag  coefficient,  neglected  diffraction  loading 
altogether,  and  also  used  the  ambient  air  density  instead  of  shock  gas  density  in 


calculating  gas  velocity  behind  the  shock.  Mercx's  model  is  noteworthy  because  it  is 
commonly  used  by  the  oil  and  gas  industry,  especially  in  Europe. 

The  effects  of  the  deficiencies  in  the  two  models  are  noted  in  subsequent  discussion.  Eqs.  (1) 
through  (12),  as  formulated,  resolve  such  deficiencies  and  constitute  the  Improved 
Displacement  Velocity  (IDV)  model.  Numerical  solution  of  the  IDV  model  equations  was 
accomphshed  by  means  of  the  EAI  SWEEP  software  for  pressures  ranging  from  .01  to  1000 
psi  and  impulses  ranging  from  1  to  10^  psi  ms.  The  results  were  then  plotted  in  the  form  of  a 
P-I  diagram  for  purposes  of  comparison  with  the  results  produced  by  Baker’s  and  Mercx’s 
models,  as  described  in  Section  3.0 

3.0  COMPARISON  QE  RESUETS 


In  order  to  compare  the  results  produced  by  the  three  models  a  series  of  calculations  were 
carried  out  for  the  case  of  6-foot  tail,  180-pound  man.  The  resulting  pressure-impulse 
diagrams  based  on  Baker's  model,  Mercx's  model,  and  the  IDV  model  are  presented  in 
subsections  3.1,  3.2,  and  3.3  respectively.  The  differences  in  the  predicted  vulnerabihties  are 
discussed  in  subsection  3.4. 

3.1  BAKER'S  MODEE  RESUETS 

Based  on  Baker's  model,  the  four  displacement  velocities  given  in  Table  1  can  be  represented 
by  a  family  of  four  curves  on  a  pressure  impulse  diagram,  as  shown  in  Eigure  1 .  In  this 
diagram  the  abscissa  is  peak  overpressure,  Pq,  and  the  ordinate  is  positive  phase  impulse,  Ipp. 
Examination  of  Eigure  1  reveals  the  effect  of  neglecting  the  changing  displacement  velocity  as 
a  function  of  time  in  calculating  drag.  Erom  the  basic  principle  of  compressible  fluid  flow  the 
gas  velocity  can  be  expressed  as  a  function  of  overpressure  as  given  by  Eq.  (7).  Eor  typical 
atmospheric  conditions  [17]. 

Pqo  =  14.6941  psia 

Too  =  531.69  °R 

Roo  =  53.36  ft-lbf/(ft-lbm°  R) 

With  these  values  introduced  into  Eq.  (7),  the  maximum  gas  velocity  behind  the  shock 
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Figure  1.  Pressure-Impulse  Diagram  for  Human  SkuU  According  to  Baker's  Model  [7] 
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can  be  represented  as  fiinetion  of  the  peak  overpressure,  Pq,  as  shown  in  Figure  2.  Now 
the  maximum  displaeement  veloeity  eannot  exeeed  the  maximum  gas  veloeity  behind  shoek, 
espeeiaUy  for  low  pressures  and  large  impulses.  Thus,  based  on  Figure  2,  for  eaeh  of  the  four 
displaeement  veloeities  there  is  a  minimum  overpressure  below  whieh  the  dis-plaeement  veloeity 
eannot  oeeur.  These  minimum  overpressures  are  tabulated  in  Table  2. 


As  indieated  in  Figure  1,  when  these  minimum  overpressures  are  taken  into  aeeount,  the  upper 
one-third  of  Baker's  eurves  are  physieaUy  impossible.  Thus,  Baker's  curves  greatly  over¬ 
estimate  the  injury  to  the  skull  for  low  pressures  (<  .5  psi)  and  large  impulses  (>  20,000 
psi-ms). 


Table  2.  Minimum  Peak  Overpressures 


Displaeement 
Veloeity  (ft/see) 

10 

13 

18 

23 


Minimum  Peak 
Overpressure  (psi) 

.1829 

.2382 

.3307 

.4237 


3.2  MERCX'S  MODEL  RESULTS 


In  the  case  of  Mercx,  the  four  displacement  velocities  can  also  be  represented  as  a  family  of 
four  curves  on  a  pressure-impulse  diagram,  as  shown  in  Eigure  3.  Comparison  of  these  curves 
with  Baker's  curves  reveals  a  general  difference  in  shape.  Eor  equal  displacement 
velocities  close  inspection  reveals  that  the  curves  do  not  intersect,  with  Mercx's  curves  always 
lying  above  and  to  the  right  of  Baker's.  No  physically  obvious  discrepancy  can  be  detected  by 
simple  inspection,  but  the  fact  that  Mercx's  curves  share  no  common  points  with  Baker's  is 
interesting.  Eor  low  pressures  the  difference  in  the  curves  results  from  the  deficiency  in 
Baker's  drag  calculations.  At  higher  pressures  the  difference  is  caused  primarily  by  Mercx's 
use  of  an  unreahstically  low  drag  coefficient  (-.75)  and  his  neglect  of  diffraction  loading  on  the 
body.  For  these  higher  pressures  Mercx's  curves  tend  to  underestimate  the  vulnerability 
of  the  skull. 

3.3  IDV  MODEL  RESULTS 

By  means  of  the  IDV  model,  the  four  displacement  velocities  produce  four  curves  in  the  P-I 
diagram  shown  in  Eigure  4.  Close  inspection  reveals  that  a  low  pressure  and  larger  impulse  the 
IDV  curves  closely  match  Mercx's  while  at  higher  pressure  sand  low  impulse  the  IDV  curves 
agree  fairly  well  with  Baker's.  This  pattern  is  clearly  shown  in  Eigure  5  for  the  18  ft/sec 
displacement  velocity  curve.  Thus,  the  IDV  curves  at  low  pressures 
generally  result  in  lower  vulnerability,  when  compared  with  Baker's,  and  at  high 
pressures  result  in  greater  vulnerability,  when  compared  with  Mercx's.  These  results  are 
completely  consistent  with  the  more  rigorous  formulation  of  the  problem  associated  with  the 
IDV  model. 

In  order  to  increase  their  utihty  the  IDV  curves  for  the  skuU  have  been  incorporated  into  the 
Vulnerabihty  Assessment  of  Structurally  Damaging  Impulses  and  Pressures,  (VASDIP)  3.0 
software  [18].  Similar  curves  for  eighteen  other  components  of  the  human  body  are  also 
included  in  the  software. 

3.4  DIEEERENCES  IN  VULNERABILITY  PREDICTIONS 

To  demonstrate  the  vulnerabihty  differences  produced  by  the  three  models,  a  series  of 
calculations  were  carried  out  along  the  IDV  18  ft/sec  curve,  corresponding  to  50%  injury. 
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Figure  4.  Pressure-Impulse  Diagram  for  Human  SkuU  According  to  IDV  Model 
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Figure  5.  Comparison  of  Pressure-Impulse  Diagram  for  18  ft/sec  Displacement  Velocity 


The  results  are  shown  in  Table  3  and  Figure  6.  Again  these  results  reinforce  the  statements 
already  made  concerning  the  inaccuracies  associated  with  both  Baker's  and  Mercx's  models. 
Agreement  between  all  three  curves  is  not  achieved  at  any  point  along  the  plot.  While  these 
differences  do  not  appear  to  be  that  great  when  comparing  the  famihes  of  curves  on  P-I 
diagrams,  the  differences  in  vulnerability  as  given  in  Table  3  and  Figure  6  differ  by  as 
much  as  80%. 

Table  3.  Comparison  of  Predicted  Vulnerabihty  Along  IDV  18  ft/sec  Curve 


Overpressure 

Impulse 

(psi) _ 

(psi  msl 

0.5 

83,619 

1 

11,445 

2 

3,758 

5 

1,162 

10 

511.2 

20 

222.1 

50 

72.53 

100 

39.59 

200 

26.48 

500 

20.07 

1000 

17.37 

Vulnerabihty  (%) 


Baker 

Mercx 

IDV 

100 

50 

50 

100 

19 

50 

66 

18 

50 

31 

7 

50 

36 

0.7 

50 

— 

0 

50 

— 

0 

50 

— 

0 

50 

— 

— 

50 

— 

— 

50 

_ 

_ 

50 

4.0  CONCLUSIONS 


The  two  most  commonly  used  models  for  skuU  vulnerabihty  [7,  10]  each  suffer  from  significant 
deficiencies  which  result  in  vulnerabihty  values  which  consistently  differ  from  one  another  from 
25  to  80%.  The  more  recent  IDV  model  was  designed  to  overcome  the  deficiencies  associated 
with  the  two  earher  models,  and  the  results  demonstrate  that  the  IDV  model  serves  as  a  bridge 
between  the  earher  two.  Because  of  the  importance  of  skull  injuries  in  predicting  blast 
effects  on  the  human  body,  the  use  of  the  IDV  model,  as  incorporated  into  the  VASDIP 
3.0  software  [18],  should  significantly  improve  the  accuracy  of  such  injury  predictions. 
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Figure  6.  Predicted  Vulnerability  Along  IDV  18  ft/sec  Curve 
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